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We investigate spontaneous four wave mixing (SFWM) in a single-channel side-coupled integrated 
spaced sequence of resonators (SCISSOR). Analytic expressions for the number of photon pairs 
generated, as well as the biphoton wave function (joint spectral amplitude) describing the pairs, 
are derived and numerically computed for different pump pulse durations and numbers of ring 
resonators. In the limit of a long input pump pulse, we show a strong analogy between super- 
linear scaling of generation efficiency with respect to the number of rings in the structure and Dicke 
superradiance. More generally, we discuss in detail the factors that influence the shape of the 
biphoton wave function, as well as the conditions for observing super-SFWM. 


I. INTRODUCTION 

The last few years have witnessed a burgeoning interest 
in integrated quantum nonlinear optics [lH8|. Nonlinear 
optical processes, such as spontaneous parametric down 
conversion (SPDC) and spontaneous four-wave-mixing 
(SFWM), have been explored as reliable and efficient 
sources of nonclassical states of light in integrated semi¬ 
conductor structures. The development of such sources 
is one of the major challenges in implementing a full inte¬ 
gration of quantum optics on-a-chip Q, which holds the 
promise of a host of applications ranging from quantum 
information processing to quantum computing. There 
is also an interest in integrated quantum optics beyond 
developing technological applications, as quantum optics 
has always been an important resource for the study of 
fundamental concepts in quantum mechanics, including 
entanglement and weak measurement m- More re¬ 
cently, integrated optics has been used to demonstrate 
fundamental results such as Anderson localization in the 
quantum regime [l^ and boson sampling [l3l - [l5j |. 

Both SPDC and SFWM have corresponding classi¬ 
cal processes, respectively difference frequency genera¬ 
tion and stimulated four-wave-mixing. Most of the con¬ 
ditions that lead to high efficiency for these classical 
phenomena are those required to achieve l^h efficiency 
for the corresponding quantum processes |16l - [lq . Thus, 
the design of integrated devices to enhance the nonlin¬ 
ear light-matter interaction in the quantum regime can 
take inspiration from recent results concern ing classical 
frequency conversion in integrated devices |19l - l^ . and 
from the know-how developed over more than fifty years 
of classical nonlinear optics. Similarly, it makes sense to 
investigate parametric fluorescence in devices that have 
been initially proposed for applications in linear optics, 
as many of them have peculiar properties that could be 
exploited in the quantum nonlinear regime. In part, this 
was the motivation for our recent work on SFWM in 
single-channel side-coupled integrated spaced sequences 


of resonators (SCISSORs) [1^. Known for years as linear 
components in integrated optics [1^ , we showed that 

these systems have interesting properties in the nonlin¬ 
ear quantum optics regime: for continuous wave (CW) 
pump excitation, the probability of photon pair produc¬ 
tion by SFWM can scale quadratically with the number 
of resonators, a phenomenon we named super-SFWM. 

Superradiance, or the phenomenon of several emitters 
acting in a coherent way that leads to an enhanced emis¬ 
sion rate, was first explored by Dicke in his classic 
paper “Coherence in Spontaneous Radiation Processes,” 
in relation to the spontaneous emission of a gas of atoms. 
In the present work we cast superradiance in a more gen¬ 
eral light, and consider super-SFWM as an instance of 
quantum nonlinear optical processes acting coherently. 
While such coherence is quite general, and can be ob¬ 
served for SFWM in a variety of structures, here we use 
the simple example of a SCISSOR, which we believe best 
illustrates the connection between the superradiance of 
super-SFWM and that of Dicke’s work. Extending our 
earlier work on generation efficiency with CW pump¬ 
ing (^ . here we formalize the aforementioned connec¬ 
tion, relax the CW approximation, and explore the spec¬ 
tral properties of the state of generated photons in addi¬ 
tion to its generation efficiency. 

We begin in Section |TI] with a simple calculation based 
on Fermi’s Golden Rule. While this is strictly limited to 
identifying the rate at which pairs are generated, it has 
the advantage of explicitly showing a parallel between 
coherence in parametric fluorescence and coherence in 
spontaneous emission. We start with a general discussion 
of the simple theory of superradiance, and identify the 
parallels between Dicke’s example of atoms in a gas and 
the example of rings in a SCISSOR structure. This sim¬ 
ple calculation shows that super-SFWM can arise from 
the constructive interference of the probability amplitude 
of generating a photon pair in each ring. In Section IIIII 
we show how, in the limit of a long pump pulse, general 
expressions for SFWM in a SCISSOR reduce to the sim- 
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pier Fermi’s Golden Rule result of Section im details of 
the application of the asymptotic fields formalism [2^ to 
a SCISSOR, the long input pump pulse limit applied to 
these expressions, manipulations involved in the connec¬ 
tion with Dicke superradiance, and the derivation of the 
maximum number of rings that are expected to behave 
coherently are relegated to Appendices 0 [B] [Cl and iDl 
respectively. 

Armed with a general description of SFWM in a SCIS¬ 
SOR, even in the presence of pump pulses with duration 
comparable or shorter than the dwelling time of a pho¬ 
ton in each ring resonator, we go beyond the long pulse 
limit in Section ITVl and study the time-energy correlation 
of the generated photon-pairs in general. We evaluate 
the biphoton wave function (joint spectral amplitude) 
that characterizes them, revealing how the interference 
depends both on the spectral features of the pump pulse, 
and on the geometrical parameters of the structure. In 
particular, this suggests a strategy for the design of ad 
hoc biphoton wave functions, in which the quantum cor¬ 
relations of the generated photon pairs are tailored by ad¬ 
justing the interference between different sources driven 
by the same pump pulse. We examine how the gener¬ 
ation rate of photon pairs depends on the properties of 
the pump pulse, and particularly how the number of rings 
that radiate coherently depends on that pump pulse. We 
also present some numerical results for a specific struc¬ 
ture, and illustrate the biphoton wave functions that re¬ 
sult as the number of rings and the pump pulse duration 
are both varied. This allows us to identify the conditions 
that must be satisfied to observe super-SFWM in a given 
structure, and how they depend on geometrical and ma¬ 
terial parameters. Some considerations, and connections 
to other third-order nonlinear optical processes, are dis¬ 
cussed in Section |Vl and our conclusions are presented in 
Section |Vl| 


II. SUPERRADIANCE: ATOMS AND RINGS 


To ground our analysis of super-SFWM in something 
familiar, we begin with a simple Fermi’s Golden Rule 
calculation that establishes a connection between super- 
SFWM and Dicke superradiance. Discussions of super¬ 
radiance typically involve an output, or “radiation” sub¬ 
system with Hilbert space "Hr, consisting of at least part 
of the radiation field, and a “pump” subsystem, consist¬ 
ing of N identical or nearly identical components, with 
Hilbert space 'H-p. The full Hilbert space is a direct prod¬ 
uct of these two subsystem Hilbert spaces, 

■H = nn^np. ( 1 ) 

In Dicke’s work [2^ the output subsystem is the full ra¬ 
diation field and its initial state is the vacuum, while the 
pump subsystem is a collection of atoms and the initial 
state is an energy eigenstate. We slightly generalize the 
usual elementary treatment of superradiance by consid¬ 


ering an initial state of the form 

|iA(0)) = |vac) 0 l^-p), (2) 

where |vac) indicates the vacuum state of the radiation 
subsystem, and I'I'p) is now an arbitrary initial state of 
the pump subsystem; for the moment we do not consider 
how the components of the pump subsystem were initial¬ 
ized in the state I'I'p). The Hamiltonian for the system 
can be written very generally as 

H = Hp + Hp + V, (3) 

where Hp and Hp are respectively the Hamiltonians of 
the isolated radiation and pump systems, and V describes 
the interaction between them. We will see that in our 
examples that the salient part of the interaction V can 
be written as 

V = -J2cra<inlVm+B^-C., (4) 

m,q 


where the Cmq are constants, is a raising operator 
acting on kets in the radiation Hilbert space, Vm is a 
destruction operator acting on kets in the pump Hilbert 
space, and H.c. denotes Hermitian conjugate. Here m la¬ 
bels the component of the pump system, m = 1, 2,..., A^, 
and g is a parameter characterizing the output system. 
For an initial state ([2]) the ket will evolve to |'i/’(i)) at time 
t, and we seek the probability ^r(t) that an appropriate 
measurement performed on the radiation system would 
find it in an energy eigenstate \ijjr) of Hp, independent 
of the state of the pump system, 

Vr{t) =^\{{tj;p\iS: {A\)\ll^{t))\^ , (5) 

P 


where the |'0p) are energy eigenstates of Hp. Applying 
perturbation theory to calculate |V’(^))i a calculation un¬ 
der the usual assumptions of Fermi’s Golden Rule leads 
to the result 


^r{t) 

t 


=—6(hi0r — hn) 
h 

X ^ (4'pIQ'^ (m,r) Q(m',r) |4'p) , 

m,m' 


( 6 ) 


where Hplijjr) = hurlilJr), we take the energy of the 
vacuum state of Hp to vanish, and we have assumed 
^iHpt/h'p^^-iHpt/h _ some frequency D; 

this can easily be generalized, but it will suffice for the 
examples of interest here. The operators Q (m, r) = 
Cmq {i^r\Hy[vsic)'Prn are Operators over the pump 
Hilbert space that depend parametrically on the final 
state |'0r) of the radiation system being considered; it 
is understood that is to be summed over all possi¬ 
ble final states IV'r), with that sum then passing to an 
integral in the continuum limit. Finally, we write (|6]) as 


^rit) 

t 



hQ) 


l^-p) 


2 


( 7 ) 
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and see that the initial radiation rate depends on how 
lEm 2 scales with increasing N. If it scales 

as -/V - that is, there is constructive interference between 
all the kets Q {m, r) I'I'p) for different m - then the initial 
radiation rate scales as . 


A. Atoms 

We first recall how these results apply to the usual 
problem of Dicke superradiance, where the pump system 
consists of N identical two-level atoms and the output 
system consists of modes of the radiation field. Here we 
have 


m 

( 8 ) 

where bq is a destruction operator for a mode of the radi¬ 
ation field, and we take q to label both the wave vector k 
and the polarization type; huj is the transition energy of 
the atoms, and = (|e)„ (e|„ - |5)™(<?L)/2, where 
|c)m ^cid \g)j^ are respectively the excited and ground 
states of atom m. The interaction Hamiltonian in the 
dipole approximation is given by 

^ = -E®™-E(rJ, (9) 


where 2)™ = iZmC™ -I- H.c. is the dipole moment op¬ 
erator of the atom and E(r^) is the electric field 
operator at the position of the atom. Here Um is 
the vector electric dipole matrix element for the transi¬ 
tion, cr™ = \g)^ (e|^ is the lowering operator for the m*** 
atom, and E(r) = iSqeqbqe^^’’^ + H.c. This operator is 
defined by the polarization of a plane wave mode and 
the electric field per photon £q = yjhckq/2eoV, where 
V is the normalization volume. For the energy conserv¬ 
ing part of the interaction we find Vm = o™, 
and Cmq = . To calculate the probabil¬ 

ity per unit time that a photon is emitted into the state 
\ipr) = bl |vac) = |r), we identify H = a), and find 


Vrjt) 

t 



tUjj) 


'^Q{m,r) l^-p) 


2 


( 10 ) 


from with 


g(m,r) = (11) 

For Dicke’s example [1^ of N identical atoms localized 
within a wavelength of light, the Cmr are independent of 
TO, Cmr —t Cr, and the sum over m of the Q {m,r) is pro¬ 
portional to the sum over m of the cr™. The states |4'p) 
that lead to superradiance are the states \JM), which 
are completely symmetric with respect to exchange of 


the atoms, and where J + M atoms are in the excited 
state, and J — M atoms in the ground state. For these 
states we find 

= ‘^-d{njJJr-tluj)\Cr\^ {J + M){J-M+l). (12) 

This takes its largest value when J = N/2 and M = 0, 
corresponding to half of the atoms in their excited states, 
and scales there as N'^. 

0T) O O 

-^^-//-^^-► 

FIG. 1. The SCISSOR structure. All N rings have radius R 
and are equally spaced from their neighbors by a distance A. 
We position the first ring at zi , the second at Z 2 = zi + A, 
etc.; thus zn — zi = {N — 1) A. 


B. Rings 

Next we turn to the SCISSOR structure composed of 
N rings coupled to a single bus waveguide. For simplic¬ 
ity, we assume that all of the rings are identical, with 
an equal distance A between each ring and the follow¬ 
ing ring, until the ring is reached (see Fig. [1}. To 
look at super-SFWM in this configuration we take 'Hp 
to be the Hilbert space of the radiation field at pump 
frequencies, and Hr the Hilbert space of the radiation 
field at the signal and idler frequencies. Here Hr in¬ 
volves only part of the full radiation field, but we refer to 
that part as the “radiation subsystem.” In SFWM the 
pump frequencies are typically in the neighborhood of 
one of the resonances wp of the ring resonators, and the 
signal and idler frequencies are in the neighborhood of 
resonances at higher and lower frequencies respectively, 
wi and ws (see Fig. [2]). Light at all these frequencies can 
in reality pass between the bus waveguide and the rings. 
However, for our calculation of super-SFWM in this Sec¬ 
tion we neglect any coupling of light between the rings 
and the bus waveguide at the pump frequencies. Then 
as an initial state one can imagine pump light “loaded” 
into the rings, and no light present at signal and idler 
frequencies; at the moment we do not concern ourselves 
with how this initial state of the rings was prepared, just 
as in the elementary calculation above of Dicke’s system 
we did not consider the preparation of the initial state 
of the atoms. In the SCISSOR example our initial state 
can then also be taken as where now |vac) refers to 
the vacuum state of the radiation field at signal and idler 
frequencies, and I'I'p) identifies the initial state of the 
pump light in the rings. In this idealization, were there 
no SFWM the pump light would remain in the rings for¬ 
ever, just as if there were no spontaneous emission the 
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atoms introduced by Dicke would remain in their initial 
state forever. We assume that the resonance modes of 
the rings of interest are those where light is propagating 
counterclockwise within the rings. 


physical system; we quantize the output system over the 
length L under the usual assumption of periodic bound¬ 
ary conditions, and will ultimately take C ^ oo. The 
Hamiltonian for the output light is 



Hr = ^ h^{k)h\bk + ^ fvjj{k)h\bk: (15) 

kG{ks} k^{ki} 


where 



Skk', with {fcs} and {ki} respectively 


the sets of signal and idler frequencies under the appro¬ 
priate resonances. The wave number k here labels the 
different modes, with frequencies uj{k). Since the signal 
and idler light propagates through the rings as well as in 
the bus waveguide, and since we are interested in signal 
and idler light being generated and leaving the system, 
we take these modes as “asymptotic-out” states of the 
radiation field [2^, familiar from scattering theory [s^- 
In this example, where we are interested in calculating 
signal and idler light propagating towards 2 = -l-oo, the 
asymptotic-out states are chosen so that for z greater 
than the position of any of the rings but within the cen¬ 
tral normalization length the fields are propagating with 
wave number k and normalized as would be the fields 
propagating in an isolated bus waveguide; their behavior 
for smaller 2 ; along the waveguide, and within the rings, 
then follows from the Maxwell equations. From a stan¬ 
dard introduction of these asymptotic-out fields [ 2 ^ we 
write the displacement operator associated with them as 


FIG. 2. Absolute value squared of the field enhancement fac¬ 
tor in the ring (1391) . with A the FWHM of the resonances and 
FSR the free spectral range. 


Within this framework the Hamiltonian for the pump 
system describes the resonator modes in the rings at the 
resonance where we consider the (assumed isolated) res¬ 
onators are “loaded” with pump light, 

Hp = bwp'^a'l.am, (13) 

m 


where henceforth wp is the pump resonance of interest 


and 


1 ^ 


t 


= Sr, 


Quantizing in terms of the elec¬ 


tric displacement and magnetic fields (29|, the displace¬ 
ment field operator in the rings at the pump frequency 
can then be written as 


Dp(r) = 




hijp 


H" H-C., 


(14) 


where Diso(r, m) is a mode field that exists only in ring 
m, and describes the counterclockwise propagation of 
light (see Fig. [3]). 

Through SFWM this will generate light at signal and 
idler frequencies that will ultimately leave the rings and 
propagate towards 2 ; = 00 in the bus waveguide. To 
describe this radiation, which here is the “output sys¬ 
tem,” we introduce a normalization length C of the chan¬ 
nel, much longer than the size (about (iV — 1) A) of the 


Do(r) + H.c. 

+ VZ S “‘(r) + Hx., 

( 16 ) 

(see Fig. |31) where S and I refer to signal and idler 
photons respectively, and where the subscript R on the 
asymptotic-out mode D^Y“°“*(r) indicates light leaving 
the system to the right; while such an expansion is usu¬ 
ally done where k ranges continuously (29|, the overall 
factor y^27r/£ here accounts that for the moment we are 
quantizing within a normalization length C. 

o 0 0 O.,^ 

FIG. 3. A sketch of the Dicke picture of super-SFWM in a 
SCISSOR, with with strong pump light shown confined to the 
rings and asymptotic-out signal and idler states. 

The nonlinearity responsible for SFWM is character¬ 
ized by the usual third order susceptibility tensor X 3 [SH ■ 
Writing the nonlinear Hamiltonian in terms of the dis¬ 
placement fields instead [s^, the appropriate V appear- 
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ing in ® is 

^ = -irE / drrf'(r)Z?*(r)i^^(r)7^^(r)Z?'(r), 

(17) 

where D(r) is the total displacement operator, and we 
neglect dispersion in the nonlinearity; r 3 '^^^(r) is related 
to (y), and the integral over Vm indicates an integral 
over the ring. Using m and in (II3, we keep 
terms of the form b\^bl,.amO-m + H.c., where two pump 
photons are converted to a signal and idler photon, and 
find 


27r \ wp oJsiOi 

^ hi 

m,ks ,ki 

X F^""'>{uJs,ul^)bl^bl^amam + H.c., (18) 

where ks(i) ranges over the wave numbers associated with 
signal (idler) light, with ujg^i) = uj{ks(^i)) the associated 
frequencies, and 


F(™)(ca„caO= / drrf'(r) 


^2,asy-out 
^Rks ' 


r) 


^^asy-out(j,) m)Dl^{r, m). 

(19) 


The interaction Hamiltonian (flSl) is of the form intro¬ 
duced above ([4]), where the label q is now identified with 
the pair (fcs, fci), the wave numbers of the signal and idler 
respectively, V^n = ama-m, '^\k,,ki) = ^-iid 


^m{ks ,ki ) 




£ 


4ei 


( 20 ) 


Here H = 2u;p, and the probability per unit time that a 
signal and idler photon are emitted to form the state of 
a pair of photons IV'r) = bj. |vac), with energy huir = 
Hujs + huJi, is given by m as 


^r{t) 

t 



2huJY>) 


2 


^ Q(m,r) l^'p) 


( 21 ) 


with Q{m,r) = 

We now see how super-SFWM can arise, as well as its 
analogy with the Dicke superradiance of a collection of 
atoms. Suppose that the light at the pump frequency in 
each ring is a coherent state, with |4>p) = |{aj}), where 
{aj} denotes the set {ai, a 2 ,..., a^}, with 

Q-m |{cti}) = Q^m |{Q^j}) • (22) 


Then 


'^Q{m,r) |4'p) 



^m(ks,ki)^m 


l^p) ■ 


(23) 


While in Dicke’s example the atoms were all assumed 
to be identical and within a distance of each other that 
is much less than the wavelength of light, and thus in 
the corresponding Q (m, r) m the Cmr could be taken 
to be independent of m, here even if the rings are all 
identical the distance from one ring to any of the others 
will be much larger than the wavelength of light, simply 
because each ring itself is larger than the wavelength of 
light. Hence the will in general be different 

for each m, even though we have assumed the rings are 
identical, for the asymptotic-out field D^Y"°“*(r) in one 
ring will be different compared to that in another due 
to a phase difference arising from propagation between 
rings in the bus waveguide. Thus the Cm(ks,ki) for 
different m will in general differ by a phase. Nonetheless, 
one can imagine that if the phases of the am are adjusted 
to compensate for this, and the \am \ is the same for each 
TO, then 


^ ^m{ks ,ki)^m ^i ( 24 ) 

m 

where N is the number of rings, and super-SFWM will 
arise, with a simple description very similar to that of 
Dicke superradiance. Identifying the conditions under 
which this will occur, and constructing a more realistic 
description of the state of the pump light in the rings, is 
the goal of the following Sections. 


III. THE LONG PUMP PULSE LIMIT 


A natural limit to look for agreement between the pre¬ 
vious section and a more sophisticated treatment of non¬ 
linear quantum optics in a SCISSOR is that of a long 
pump pulse. This is because a Fermi’s Golden Rule cal¬ 
culation results in a transition rate, and as the excitation 
becomes effectively CW the probability of pair genera¬ 
tion per pump pulse can also be written as a rate. Thus 
this is where we begin a more detailed investigation of 
super-SFWM. 

In general, a SFWM process in the presence of a pump 
pulse leads to a multimode squeezed state. However, in 
the limit of a small squeezing parameter, \j3\, this state 
can be approximated as 


P f 

Itp) « I vac) + -^ J dwida;2 (j) (wi,a;2) |vac), 

(25) 

where (j)(uji,uj2) is the normalized biphoton wave func¬ 


tion (BWF), allowing interpretation of |/3|^ as the prob¬ 
ability of generating a pair of photons per pump pulse. 
Neglecting time-ordering corrections and other non¬ 
linear effects we note that this BWF provides complete 
knowledge of the full state of generated photons, and thus 
enables calculation of arbitrary field expectation values 
and correlation functions [s^. 

As we have previously sketched , and detail in Ap¬ 
pendix 1^ of this work, the normalized BWF for a SCIS- 
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SOR of N rings can be written as 
‘i'Ki'j2c?h, 


(/)(a;i,W 2 ) =- 


X 


4eo/3 
/ do;. 


a;ia ;2 


n(a;i)n(a; 2 ) 


a;(a;i + UI2 — oj) 


:</>p(w) 


J y v{ui)v{uJi + UJ2 — Lo) 

X d)p(wi + W2 — w) J(u;i, W2, w, Wl + W2 — w), 

(26) 

where ())p(u;) is the pump pulse waveform and is normal¬ 
ized according to 


J dw \(j)p (w)!^ = 1 . 


(27) 


Here J (wi, W 27 W 3 , W 4 ) is a generalized phase matching 
function defined later, the group velocities v (cUn) = 
and |q;|^ is the average number of photons 
per pump pulse. The probability |/3|^ that a pair of pho¬ 
tons is generated can be extracted from the normalization 
condition 


J |(;ii(a;i,W2)|^da;idw2 = 1, 


(28) 


in the problem, in particular compared to the dwelling 
time of a photon in a ring - the expression is well- 
approximated by 



1 97r3 


AT 2e2 



/ hwp y 

\v{ujp)J 

uj(2u!p — uj) 
v(uj)v(2ujp — uj) 


I J {uj ^ 2ujp 


w,a;p,wp)| , 

(32) 


where the integration over w now ranges only over the 
frequency range of the signal (see Fig. HD. For a fixed 
number of rings the only dependence of the right-hand- 
side on the pump pulse is through its dependence on 
(AT)“^. Thus, for fixed |a|^, the probability that a pair 
of down-converted photons is generated is inversely pro¬ 
portional to the length of the pulse. This scaling arises 
simply because hujp |q;| is approximately the energy in 
the long pump pulse; the energy density (£ (energy per 
unit length) of the incident pump pulse is 


fvjjp \a'^ 

z;(a;p) AT’ 


(33) 


and is given explicitly by 


and we can write (15^ as 


3Tr\/2a^h 


4eo 
/ dw . 


dwida ;2 


wia ;2 


v{uji)v{uj2) 


Uj{uJi UJ2 — <^) 
j{uj)v{loi + a;2 — w) 


(j)p{uj) 


X (jjp {uj\ -j- UJ 2 — UJ^J{uJ\ , UJ 2 , “t“ ^2 — ^) 


(29) 


The long pump pulse limit of these expressions allows a 
strong connection with the Fermi’s Golden Rule result of 
the previous Section. 

We first consider a pump pulse of length AT idealized 
as a “top hat function” with a center frequency wp. 




AT AT 

< t < 

2 2 


=0 otherwise. 


(30) 


We take ujp to be one of the ring resonances. The spectral 
amplitude (j)p{^) is then given by 


(j)p{uj) 


1 . / W — UJp 

= , smc —-- 

V ttAuj V Auj 


(31) 


where Auj = 2/AT and sinc(a;) = (sina;)/x, satisfying 
( 071 ) . In Appendix|B]we show that for Auj sufficiently nar¬ 
row - i.e. AT long enough compared to other time scales 


AT 


97r3 

Z£g 


uj{2ujp uj^ I T/" o ^l^ 

auj , ' —- Tj\J{uj,Zujp — uj^ujp,ujp)\ . 


v{uj)v{2ujp — uj) 


(34) 


showing that the rate of photon pair production depends 
on the square of the energy density of the pump pulse, as 
expected for spontaneous four-wave mixing. At the same 
level of approximation we find in Appendix |B] that the 
expression for the joint spectral density is 


|</(wi,u; 2 )|^ 


v{uJl)v{uj2) 

X \J {uji,uj 2 ,ujp,ujp)\^ 5{2ujp 


UJl - UJ 2 ), 

(35) 


where wi and uj 2 can range over both signal and idler 
frequencies, and A/” is a constant to guarantee the nor¬ 
malization condition (l28l) . 

The presence of the energy-conserving Dirac delta 
function in (IMll . and the prediction (|Mll of a rate of pho¬ 
ton pair emission, suggest that the long pump pulse limit 
identified here should be described by the simple Fermi’s 
Golden Rule argument given in Section III Bl However, 
to make that connection we have to relate the pair pro¬ 
duction rate (IMl) not to the square of the energy den¬ 
sity in the pump pulse, but to the square of the energy 
Swp that the pump pulse “loads” into each ring (re¬ 
call (HI,®), since it is these that appear in m, 
(1231) . We do this in Appendix o and indeed find that 
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(l34ll is equivalent to the result (1^ for the Fermi’s Golden 
Rule calculation. 

What is yet to be established is how the rate of pho¬ 
ton pair emission scales with the number of rings. We 
address that issue in the next Section, for both the long 
pump pulse limit that we have shown is described by the 
simple Fermi’s Golden Rule calculation, and for shorter 
pump pulses where it is necessary to use the more general 
formalism. 


IV. THE COHERENCE OF THE RINGS 


To evaluate the pair production rate, either in the long 
pump pulse limit (1341) or more generally (j29l) . we now 
examine the generalized phase matching function 


J(wi, a;2, W3, W 4 ) 


sin 


HN 


sinf 


X [k{uji), k{u}2), k^uJs), fc(w4)], 


(36) 


which is derived in Appendix El the phases 
X (wi, 012 , W 3 , 014 ) and ^ (wi, W 2 , W 3 , W 4 ) are discussed 
below. The two crucial components of dSH) are 
fc(w 2 ), fc(w 3 ), fc(w 4 )], the contribution from a 
single ring, and the factor involving ^(wi, 012 ,^ 3 , 014 ) 
that identifies to what extent the contributions from 
the different rings add in phase. If /i = 0 then there is 
complete constructive interference, or \ J\^ = N'^ , 

resulting in an N'^ scaling of the pair generation 
rate. We begin by simplifying the expression for 
f^^k{uji),k{uj 2 ),k{uj 3 ),k{uj 4 )], since the resonance 
structure of this factor in J(a;i, 012 , 013 , 014 ) identifies 
approximations that can be made in evaluating the other 
factors. 

As detailed in Appendix!^ we can write the contribu¬ 
tion from a single ring as 


[fc(wi), fc(w2), k{uJ3),k{uj4)] 


= I 


Dl 


(!■) (!•) 


(3) VC 

''k{uj3) ^ -^k{uj4) 

where the components of the displacement fields 


X Dt,..., (r) dL. (r), 


Dk.) 


(37) 


(38) 


are those inside the ring and C ranges from 0 to L Here 
the 


(39) 

are field enhancement factors with k and cr the usual 
cross- and self-coupling coefficients [s^, and I = 2ttR 
the circumference of the ring with radius R. We neglect 


the weak wave number dependence of the I?(r) in (1371) : 
more general expressions are given in Appendix The 
integral in (1371) is over a ring imagined just above the bus 
waveguide at z = 0 (see Fig. SI). Similarly, we neglect 
any weak dependence of a and n on fc(w) in the field 
enhancement factors (w) and consider the physically 
relevant case of weak coupling when 

|k| <|cr|, 

kl ^1- (40) 

We consider both k and cr positive, though note that 
generalizing this is straightforward. Resonances occur at 
frequencies ujm where 

k(ujM)l = 27rM, (41) 

with M an integer, and for uj close to ojm the field en¬ 
hancement factor (15^ simplifies to 

^ (‘^) ^ \J 1-a {lo-ujm)-i{A/2)^ (^^^ 


where 


2(1 - a)vg 


(43) 


is the FWHM of the resonance in |5'(w)|^ at ujm- We 
sketch |g' (w)l^ near our resonances of interest in Fig. [2] 
In arriving at this expression, we have assumed that uj 
is sufficiently close to lom that group velocity dispersion 
can be neglected, with the group velocity taken to be Vg . 
Indeed, we make the stronger assumption that over the 
frequency range including the pump, signal, and idler res¬ 
onances we can neglect group velocity dispersion, writing 
v(oj) = Vg for all such frequencies, where 


k{uj) = k{ujM)+ - —(44) 
'^g 

with UJM any resonant frequency in this range. This ap¬ 
proximation is justified for realistic structures, as detailed 
below [ 2 ^. With this approximation the free spectral 
range is 


^ cum-h - 

27r 



(45) 


In the weak coupling limit (l40l) it is also true that 


^ « FSR (46) 

ZTT 


(see Fig. [2|). Taking a signal ring resonance frequency 
uis and an idler ring resonance frequency uji of interest, 
with frequencies smaller and larger than the pump ring 
resonance frequency wp, respectively, such that 2 fc(a;p) — 
fc(a;s) — k{uJi) = 0 , (HI)) tells us that 


2 ajp — ws ~ wi = 0, 


(47) 












for the resonance frequencies themselves. For wi close to 
uis, W2 close to w/, and both wa and UJ 4 close to wp, (EZl) 
reduces to 


neous four-wave mixing will indeed be superradiant with 
a generation rate of photon pairs scaling as N'^; the final 
result to be used in (IMl) is 


[fc(wi), fc(w2), fc(w3), fc(w4)] 

_ / 2 y A/2 A/2 

VI - crj (wi - ws) - i(A/2) (a;2 - oJi) - i(A/2) 
A/2 A/2 

(w3 — wp) — i(A/2) (tiJ4 — wp) — i(A/2) 

X I j^r^^{r)[V\r)V^ir)]*V\r)V\r) 

X ("43^ 

We use this expression below to evaluate J(wi, W2, W3, ^4) 
for various pump pulse durations. 



FIG. 4. A single ring of radius R positioned at « = 0 with 
transverse field profiles in the ring and channel indicated by 
Gaussians. 


I J(a;, 2wp — (jJjWp, wp)/ 



(A/2)^ 

ccs)2 + (A/2)2 


2 


X 


/ 


Hr ^ 

(r) [V\r)V^ (r)] * (r) , 

(49) 


and the joint spectral density then follows from (1551) . Of 
course, for signal and idler sufficiently far in frequency 
from the pump, group velocity dispersion need be taken 
into account and the full coherence of the spontaneous 
emission rate will be lost, but that can be described by 
a straight-forward generalization of the above equations. 
We have considered the effect of group velocity dispersion 
earlier , and found that for a typical dispersion in sil¬ 
icon structures, ^ ~ 0.01 ps^/m [—8 ps/(nm-km)], about 
a thousand rings would be necessary for degradation in 
the scaling as N'^ to set in. Thus we are indeed justified 
in neglecting group velocity dispersion throughout this 
Section. 

More interesting are the generation rate and joint spec¬ 
tral density that will result if we consider shorter pump 
pulses. For that we return to the more general equations 
(1551) and (1551) . for and |/3|^, respectively. In 

these expressions uii and UJ 2 can both range over signal 
and idler photon frequencies, but we focus on the contri¬ 
butions where the first of these is close to a signal ring 
resonance frequency ws and the second is close to an idler 
ring resonance frequency wi; the term where uii and W2 
are interchanged can be considered similarly. Then the 
quantity of interest is J(oji,uj 2 ,uj 3 ,lo 4 ) where L 03 and 0 J 4 
are pump pulse frequencies, and the energy conservation 
constraint implicit in the general equations (|56l) and (l29l) 
identifies 


Within the simple Fermi’s Golden Rule description 
of Section El the interference between the biphoton 
amplitudes generated for each ring depends on the 
Cm{ks,ki)a‘^ of (1541) . and whether or not they add con¬ 
structively depends on the relative phase of the different 
amplitudes am describing the excitation at the pump fre¬ 
quency in the rings. In the long pump pulse limit, the 
calculation of Section uni identifies how the rings are ac¬ 
tually “loaded” with pump light, and using (IMl) we can 
examine what kind of radiance is actually exhibited for 
such loading. The required form of the generalized phase 
matching function J(a;, 2a;p —w, u;p,a;p) can be extracted 
as a special case of the more general problem considered 
below, but we show in Appendix ICl explicitly that under 
our assumptions of wave number independent coupling 
coefficients and weak coupling (l40l) , at the frequencies re¬ 
quired /i in fact vanishes and thus J{u}, 2u!p — a;,a;p,a;p) 
is proportional to N. So a long enough pump pulse will 
indeed “load” the rings in such a way that the sponta¬ 


LO4 = UJl + ljJ2 — (50) 

The quantity J(a;i, a;2, a;3, a;4), with UJ 4 restricted by ((50l) 
is a richer function than the function of a single variable 
J(a;,2a;p — w,wp,a;p) required in the long pump pulse 
limit. 

Rather than considering just the “top hat function” 
pulse of Section El we consider a more general pulse 
characterized by frequency amplitudes (/)p(a;), with sig¬ 
nificant amplitudes only for w within a range d A of 
a pump ring resonance frequency ujp , 

|w — wpj < (5/2, (51) 

and phase variation small over that range. This allows 
us to consider a pump pulse short enough to explore how 
/i becomes nonzero as the pump pulse duration departs 
from the long pulse limit, and yet still long enough to take 
full advantage of the field enhancement in the ring so that 
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there is a significantly enhanced spontaneous emission 
rate. Defining 


u= (uji — ujs) + (w2 — wi). (52) 

we use the energy conservation condition (1501) and the 
resonance frequency condition (H71) to write u = UJ 3 +UJ 4 — 
2 ujp and then, noting the constraints dSU that hold for 
both pump pulse waveforms appearing in (100)) and (1001) . 
we calculate that significant contributions occur only for 
wi and W2 satisfying |m| ^ S. 

In the limit (l44l) of negligible group velocity dispersion, 
as well as that of weak coupling (IdO)) . we calculate (see 
Appendix[D| that the fi of (1501) has a typical size of about 

|/i| « (53) 

Were /i truly constant, then as the number of rings N 
were increased we would begin to see the phase-matching 
function ((501) decrease when N = Ncoh, with 

Ncoh ImI = TT, (54) 


or 


Ncoh 


TT A 

2J 


71" Tpump 
2 rdwell 


(55) 


where Tpump = 1/^ is approximately the pump pulse du¬ 
ration, and Tdweii = 1/A is approximately the dwelling 
time of light in a ring. Of course fj, is not constant, and for 
N very large we expect the probability of the generation 
of a photon pair to scale with N, as coherence between 
the rings vanishes. Nonetheless, we can take N « iVcoh 
to be a somewhat heuristic identification of when the de¬ 
pendence of that probability changes from N'^ at small 
N to at large N. 

We close this Section with some calculations of the pair 
generation rate and joint spectral densities for a typi¬ 
cal silicon-on-insulator SCISSOR structure with param¬ 
eters corresponding to operation in the telecommunica¬ 
tion wavelength range [^. We consider a structure of 
N identical rings with radius i? = 5 /rm, separated by 
a distance A = 15 /rm. We take an effective (phase) in¬ 
dex h = 2.5 and a group velocity Vg = 0.75 x 10® m/s, 
corresponding to a group index Ug = 4; we neglect group 
velocity dispersion in our calculations. We take the pump 
resonance to be the 50*^'', which corresponds to a vacuum 
wavelength of about 1570 nm; the signal and idler res¬ 
onances are taken to be those immediately above and 
below. The self-coupling a between the channel and the 
rings is given by 1 — cr = 0.0126, which leads to a FWHM 
(155)) of A = 6 X 10^° rad/s; for the pump resonance this 
gives a Q factor of Q = wp/A « 20,000. 

In the absence of group velocity dispersion and dis¬ 
persion in the nonlinear response, the relation between 
ryfc^(r) and the usual is @ 

= (») 

e^nS(r) 


where n(r) is the position dependent index of refraction. 
The standard effective nonlinear coefficient 7 that char¬ 
acterizes the third order nonlinearity is then 


7 = 


Swp 
4e3 


o'^g 


dxdy 


limn 

Xi 


[d* (x, y)d^ (x, y)] * d™ (x, y)d" (x, y) 


^(.x,y) 


(57) 


where we assume that the transverse mode profiles in the 
ring are the same as in the channel (D* (r) = d* (x, y); here 
we take n(x, y) = h in the waveguide. See Kelt et al. [l^ 
for details. We take a value of 7 = 200 (Wm)“^, in agree¬ 
ment with recent experiments on ring resonators such as 
these @ . For our pump pulses we take bandwidth-limited 
Gaussians characterized by an intensity FWHM of Tpump- 

We look first at the prediction of the generation effi¬ 
ciency, which can be characterized by the ratio |/3|^ / |q;| . 
We plot this in Fig. [5] as a function of N for different 
pump pulses. For the longest pump pulse duration con¬ 
sidered, Tpump = 1 ns, we find Ncoh ~ 100 and we would 
expect the superradiant behavior in generation effi¬ 
ciency to exist for the range of N considered in Fig. [S] 
indeed it does. For Tpump = 0.1 ns we have Acoh ~ 10, 
leading us to expect that A Ri 10 would roughly identify 
the change of scaling from A^ to A, and we see from 
Fig. [5] that it does. For Tpump = 0.01 ns the assumption 
S A used in the derivation of (1551) formally breaks 
down, but still the value of Ncoh ~ 1 would lead us to ex¬ 
pect less than superradiant behavior for all A, and that 
is confirmed by the results plotted in Fig. [S] for A more 
than a very few the scaling is linear. Finally, we note that 
while for small A generation efficiencies can be higher for 
shorter pump pulse durations than for durations in the 
long pulse limit, they scale differently with increasing N. 
Thus, although in our example the generation efficiency 
in a single ring is smallest for the 1 ns pump pulse, due 
to the quadratic scaling associated with superradiance 
this pulse will lead to the largest generation efficiency for 
larger N. 

While characterizing all of the properties of the gener¬ 
ated photon pairs requires the BWF (joint spectral am¬ 
plitude) ^(wi,a;2) (1261) . even the joint spectral density 
(JSD) |(/)(a;i, ^2)!^, where |()>(a;i, a;2)|^ da;idw2 identifies 
the probability of generating “photon 1” within dwi of 
frequency uji and “photon 2” within da;2 of frequency 
0 J 2 , contains much information about the quantum cor¬ 
relations of the generated photons. In Fig. [5] we plot a 
JSD for a pair generated by SFWM in a single ring res¬ 
onator with the parameters we have adopted, and sub¬ 
jected to a pump pulse with Tp^mp = 0.1 ns, normalizing 
the frequencies to the line width A of the resonator. The 
shape of the JSD is essentially determined by the spectral 
broadness of the pump and the resonance line width. In 
particular, since the pump pulse duration is longer than 
the dwelling time Td-^eii = 0.017 ns in the ring, along the 
direction (w2 — wi) = (wi — ws) the FWHM of the JSD 
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1 10 


N 

FIG. 5. Generation efficiency |/3|^ / \a\* vs. ring number. 

(FWHMi) is dictated by the frequency width l/rpump 
of pump pulse through energy conservation, while along 
(ijj 2 — wi) = — (wi — ujs) the JSD width (FWHM2) is de¬ 
termined by the resonance line width A, which sets the 
frequency range in which light can be efficiently gener¬ 
ated within the ring [s^. 



-1 0 1 

(cOj- C0g)/A 


FIG. 6. The JSD of a single ring for Tpump = 0.1 ns. 

It is interesting to consider the JSD for the different 
scaling regimes identified in Fig. [5] We consider the ex¬ 
amples of iV = 1,3,5. For the pulse durations adopted in 
Fig. [51 the results are shown in Fig. [71 where we observe 
that the JSD is in general a function of JV and Tpump- 
However, as would be expected, in the case of a long 


pump pulse (Tpump = 1 ns) (Fig. |7la,,d, and g) the JSD is 
essentially independent of N; here the structure is fully 
superradiant for these JV, and the contribution to the 
BWF from each ring is essentially identical. There are 
strong energy correlations due to the narrow spectrum 
of the pump pulse centered at wp, leading to a strong 
localization of the JSD about the axis oji + uj 2 = Swp, 
but due to complete constructive interference the con¬ 
tributions to the BWF from each ring are identical and 
thus no modification of the BWF occurs as more rings 
are added. 

In the case of a short pump pulse (Tpump = 0.01 ns) 
(Fig. |71:,f,i) we find a significant change in the JSD as N 
increases. For a single ring (Fig.j^:), the JSD has a larger 
FWHMi than seen for long pump pulses; because of the 
large spectral width of the pump pulse here, the FWHMi 
is limited not by that spectral width, but by the reso¬ 
nance line width A of the ring. Yet as we consider larger 
and larger N (Fig. [7]; and i), the JSD narrows, indicating 
stronger energy correlations. This can be understood as 
due to interference from the contribution of the different 
rings to the BWF. Along the axis uji + uj 2 = 2ujp we 
have u = 0 from (H71) and (1571) , and along this line where 
the energy correlations are strongest we will have ^ = 0 
from (lD3]).(lD8l). and (IDllI) . and the contributions from 
the different rings to the full BWF will add in phase. 
But moving away from this axis we will have u ^ 0, and 
thus a nonzero ^ leading to partial destructive interfer¬ 
ence from the different rings (see (ID2I) 1. Thus far from 
the superradiant regime, as we are here, as N increases 
there is a significant change in the JSD. 

Finally, for an intermediate pump pulse (Tpump = 
0.1 ns) (Fig. |7}3,e,h) we observe no significant change in 
the JSD for the small values of N chosen, as the pump 
pulse duration is sufficiently long that even at A = 5 
we are only just starting to see a decrease in scaling 
with and with it a concomitant slight modification 
in the BWF, again characterized by stronger energy cor¬ 
relations, as the number of rings increases. The scaling 
of the FWHM of the JSD along the direction identified 
by uJi — UJS = UJ 2 — uJi (FWHMi) is shown in Fig. [8] for a 
larger range of A than used in Fig. [T] 


V. DISCUSSION 

We have made a number of idealizations in our calcu¬ 
lations here. One is that the rings are identical; another 
is that we have neglected scattering losses in the rings. 
We plan to turn to including these generalizations, and 
other extensions of this work, in later publications. We 
also want to mention that such “super” behavior is not 
really that uncommon in nonlinear optics. Consider for 
example the second-order classical process of difference 
frequency generation (DFG) in a channel waveguide. For 
an injected pump power Pp at wp and an injected signal 
power Ps at ws, an idler power Pi at uJi = ujp — ojs is 
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Pulse Duration 



-10 1 -10 1 -10 1 
(co^- 03g)/A (o3^- COg)/A (co^- COg)/A 


FIG. 7. JSDs of A'^-ring SCISSORs for various pump pulse durations Tpump- From top to bottom, N takes the values 1, 3, and 
5, while from left to right Tpump takes the values 1 ns, 0.1 ns, and 0.01 ns. 


generated that is given by 
PpL^ 


Pi =Ps- 


PA 


X sinc^ < [ks{ujs) + ki{uip - ws) - fcp(wp)] 


L 


(58) 


where fcp3j(u;) are the wave vector dependences of the 
pump, signal and idler beams, A is an effective area of 
the waveguide, and P is a quantity with units of power 
that characterizes the strength of the X{ 2 ) interaction re¬ 
sponsible for the DFG [l^. Note that for perfect phase¬ 


matching {ks{uJs) + ki{uip — uJs) — kp{ujp) = 0) this quan¬ 
tity scales as signaling the presence of the coherence 
of the interaction responsible for the generation over the 
whole length of the sample, and analogous to the de¬ 
pendence in Dicke superradiance or super-SFWM. Even 
if we consider the spontaneous version of this process, 
spontaneous parametric down conversion (SPDC), the 
power generated at idler frequencies is given by 


h{ujp/2) PpL^ 

r 


(59) 


where the leading term can be thought of as the vacuum 
fluctuations in the signal power leading to the SPDC [l^ . 
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N 


FIG. 8. FWHMi of the JSD as a function of N. 

Here T is a characteristic time that identifies that power, 
and depends on the phase-matching condition at differ¬ 
ent possible idler and signal frequencies; although T can 
depend on L, that dependence is typically less than linear 
|1^ . and so scaling greater than L with channel length 
can be seen even in this simple structure. Similar super- 
linear scaling is seen in channel waveguides for four-wave¬ 
mixing processes, both stimulated and spontaneous. 

Nonetheless, we believe the SCISSOR structure we 
have analyzed here is of particular interest because it 
points the way to coherently combining contributions to 
BWFs in integrated structures. One can imagine design¬ 
ing the elements of a structure, and the way in which they 
are combined, to lead to a greater flexibility than has 
previously been achieved in determining both the BWFs 
characterizing the photon pairs generated, as well as the 
generation efficiency with which they are created. The 
SCISSOR structure is only a first example, and the sim¬ 
plest of such integrated structures that could be designed 
and fabricated to lead to new possibilities for integrated 
sources of quantum correlated light. 


VI. CONCLUSIONS 

In summary, we have investigated how the coherence 
of multiple quantum nonlinear optical sources pumped 
by the same field can lead to an alternate form of su¬ 
perradiance than that originally envisioned in a gas of 
atoms by Dicke . Focusing on the particular example 
of a SCISSOR capable of producing pairs of photons via 
SFWM, we demonstrate a direct analogy with Dicke’s 
work even though the individual microrings of the SCIS¬ 
SOR and spacings between them are both much larger 
than the wavelengths of light involved. A sufficiently 
long exciting pump pulse effectively “loads” pump fields 
in the N rings of the SCISSOR with the correct phase so 
that the spontaneous generation of signal and idler pho¬ 


tons from the different rings adds constructively, leading 
to super-SFWM - a quadratic dependence of the gen¬ 
eration efficiency on N. We have shown how this co¬ 
herence follows from a calculation using Fermi’s Golden 
Rule, which is completely analogous to the calculation 
for Dicke superradiance, and how one can see it in more 
detail by studying the actual loading of pump fields in 
the ring by the incident field. This latter approach al¬ 
lows us to study how the super-SFWM is modified by 
pump pulses of various duration. 

Additionally, there is a coherence number Ncoh that 
characterizes the maximum number of rings that can 
contribute to super-SFWM; beyond this the scaling de¬ 
creases and approaches a limit that is linearly propor¬ 
tional to N. The coherence number is proportional to 
the pump pulse length and inversely proportional to the 
dwelling time of light in a ring. For N < Ncoh the bipho¬ 
ton wave function characterizing the spectral features of 
signal and idler photons generated by the SCISSOR is es¬ 
sentially independent of V, since the contributions from 
each ring add completely coherently. However, as N in¬ 
creases beyond Vcoh the contributions from all rings add 
coherently only for very correlated pairs of signal and 
idler photons, and thus the biphoton wave function ex¬ 
hibits stronger energy correlations. 
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Appendix A: Formalism 

Here we identify the asymptotic-in states for the pump 
field and the asymptotic-out states for the signal and idler 
fields. Consider first a single ring (see Eig. 2]), which is 
located just above a bus waveguide, running along the z 
direction, at z = 0. An asymptotic-in field corresponding 
to light incident from the left [z = —oo) with wavenum¬ 
ber k is given by 

2 < 0 in the channel, 
inside the ring, 

T{k) — 2 > 0 in the channel, 

in the limit of point coupling |^ . Unlike in Section m 
here we take k to vary continuously over positive val¬ 
ues, and in this equation it can range over the frequency 
range of pump, signal, and idler light. The transverse 
mode profile in the channel is written as dfc(rj^) of the 
field, where rj_ = {x,y), and we assume that the pump. 
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signal, and idler modes of interest are all associated with 
the same transverse mode profile; this can easily be gen¬ 
eralized. The transmission coefficient T{k) is given [35| 
by 


guaranteeing that we have the desired dependence 
dfc(rj.)e*'=Vv^ for z > 0 in the channel. We could 
also construct this by finding first the asymptotic-in field 
from the right, and using the general relation 


T{k) 


a(k) - 
1 — (T(fc)e*^^ ’ 


(A2) 


where I = 2ttR is the circumference of the ring, with 
radius R, and we introduce self- and cross-coupling coef¬ 
ficients, (j{k) and K(fc), in the usual way; we take them 
to be real, with energy conservation guaranteed by the 
condition 


[D^r"(r)]*=D^r"(r) (A9) 

relating asymptotic-in and asymptotic-out fields [2^. 

We now turn to the structure of interest, with N rings 
spaced such that the minimum distance from the bus 
waveguide to ring m, m = 1,... ,iV, occurs at Zm (see 
Fig. [ID. Here the asymptotic-in field from the left is given 
by 


cr^(/c) -|- K^(fc) = 1. 
The field inside the ring is given by 

ing, , ^ in{k) X>fc(r) 


(A3) 


(A4) 


where C ranges from 0 to I and is the distance around the 
ring starting from z = 0 (see Fig. HD. In a simple model 
where the ring and channel have the same dimensions, 
X>fc(r) and dfc(rj_) would describe the same mode profile 
transverse to the direction of propagation in the ring or 
channel. Comparing the first and second expressions in 
we see that the first factor on the right-hand-side 
of ([All is a field enhancement factor, which identifies 
how much stronger the field is inside the ring than in 
the channel. It is useful to express it as a function of 
frequency, defining 


5^(w) 


in [k (w)] 

I — cr [fc (w)] ’ 


(A5) 


where fc(u;) identifies how the wavevector in the chan¬ 
nel depends on the frequency. When k{uj)l is an integer 
multiple of 27r the field enhancement factor peaks, as 
sketched in Fig. |2l thus this condition identifies the ring 
resonance, and it is by adjusting the pump frequency to 
be within one of these resonances, and looking at signal 
and idler fields at neighboring resonances, that enhance¬ 
ment of spontaneous four-wave mixing occurs [l^ . 

In this simple model the ring is an all-pass structure, 
and so |T(fe)| = I. With some algebra one can write 


T(fc) = e*®('=\ (A6) 


where 


Drr"(r) 


where 




dkir±)e''‘’‘ 

Df®(r;m) 


ink)] 

ink)] 


m dfc(rx)e 

N dfc(rj^)e 


ikz 


ikz 


z < Zi in the channel, 
inside ring m, 

Zm < z < Zm +1 in tbe channel, 
z > in the channel, 

(AlO) 


D™s(r;m)=e 


ikZn 


[T(fc)]™”^D 


™S(r-z™z), 


(All) 


and the function D™®(r) is again given by (IA4I) . In place 
of (lASp . we have 

D“r"‘(r) = T;;;4v»D‘r‘"(r) (A12) 

[T{k)] 

The field operator D (r) can be expanded in terms of ei¬ 
ther asymptotic-in or asymptotic-out fields; for the quan¬ 
tum states that will be of interest we can restrict our¬ 
selves to either using the fields D“^“™(r) or the fields 
(r); for the frequency range of the pump it is 
convenient to use the asymptotic-in fields, and for that 
of the signal and idler the asymptotic-out helds. For the 
functions dfc(r) properly normalized we then have 


D (r) = y dfc '“(r) -f H.c. 

+ / dfc + H.c. (A13) 

We use Ofc to indicate lowering operators for modes in 
the frequency range of the pump resonance, and bk to in¬ 
dicate lower operators for modes in the frequency ranges 
of both the signal and idler resonances, 


^ -1 / CT(fc)sinfc/ \ 

0(fc)=7r + /c/ + 2tan j . (A7) 

From DL^^““(r) it is particularly easy to construct 
the asymptotic-out field for fields exiting to the right, 
('") I since the ring is an all-pass structure 
(|T(fc)| = 1); we have simply 

(A8) 


Ofc, a 


t 

fc' 


h,bl, 


6{k-k'), 


[afc,afc/] = [6fe,6fc/] = 0, 


(AM) 


and all operators associated with pump fields commute 
with those of the signal and idler fields. We keep implicit 
the convention that integrals over k are over positive val¬ 
ues and those involving pump mode operators range over 
the frequency range of the pump ring resonance of inter¬ 
est, while those involving signal and idler mode operators 
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range over the frequency ranges of the signal and idler 
ring resonances of interest. 

Returning to the nonlinear Hamiltonian (HB , in 
terms of both asymptotic-in and -out fields the energy- 
conserving term involved in spontaneous four-wave mix¬ 
ing is 

V = - J dkidk2dk3dk4S{ki,k2,k3,k4)bl^bl^ak3aki 

-fH.c., (A15) 


where 


we have put 

fl = fJl-{{kn}) =(^3 + k4- ki - k2)A 

+ e{k3) + 9{k4)-eik4)-9{k2), (A22) 


and 


j {ki,k2,k3,k4) = / drr|^ 3 ^ (r) 






/,ring/ 
/C4 V 


(A23) 


S{ki,k2,k3,k4) 


1 4! 

4^^ 


fio; ( fc 1) Sw ( ^2 ) hw ( fe 3 )/iW ( ^4 ) 
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j{ki,k2,k3,k4), 

(A16) 


with 


is independent of m; it is to be evaluated for a single ring 
imagined at 2 = 0 above the bus waveguide. The factor 
/i contains the effect of the propagation of pump as well 
as signal and idler light from ring to ring, and the phases 
acquired by the propagation of pump, signal, and idler 
through the ring. Using (IA20I) in (IA18I) and performing 
the summation, we find 


j{ki,k2,k3,k4) = J dr r|^ 3 ^'(r) 


^i,asy-out / \ ^j,asy-out f \ 


The factor 4!/(2!2!) appears in (IA16|) because we here al¬ 
low ki and ^2 to span the ranges associated with both the 
signal and idler frequency regimes. Now assuming that 
the nonlinearity is only important in the rings, where 
the pump field will be strongest because of field enhance¬ 
ment, we can write 


j{ki,k2,k3,k4) ='^j^'^\ki,k2,k3,k4), (A18) 

m 


where 

j^^\k4,k2,k3,k4) = [T{k4)T{k2)f j dvr^^^v) 


sin^ 

j{ki,k2,k3,k4) = e*^^-^j"'®^(fci,fc2,fc3,/c4)- (A24) 

smf 

where 

X = xi.{kn}) ={k 3 + k4-ki- k 2 )zi + N [9{ki) + 9{k2)] 
+ iN-l)fi{{kn})/2 (A25) 

At this point we can write down the BWF describing the 
spontaneously created photons. It is usual to move from 
a k representation to an w representation; then the BWF 
(p(uji,uj 2 ), defined for positive wi and ui 2 , is the proba¬ 
bility amplitude for “photon 1” being generated within 
dwi of uJi, and for “photon 2” being generated within 
da;2 of UJ 2 in the limit of a very low probability that a 
pair is generated. The BWF is a symmetric function of 
its arguments, (j)(uJi,uJ 2 ) = (/>(w2,a;i), with both variables 
wi and W2 ranging over signal and idler frequencies. It is 
normalized according to 






(A19) 


where we have used (lAlOl) and (IA12I) . and the fact that 
[T“^(fc)]* = T(fc), which follows since T(k) is of unit 
norm. The different j^"^\ki,k 2 ,k 3 ,k 4 ) differ only by a 
phase factor, because the rings are physically the same. 
Using (lAlB) we have 

^2,^3, k4) =e^(.o(.kd+e{k2))N^i{k^+k2-H-k^)z^ 

xe*'^(™-i)rf(fci,fc2,fc3,fc4), 

(A20) 


where recall we have assumed the rings are equally spaced 
with a distance A between them. 


Zm- zi= A{m - 1), (A21) 


y |(/) (wi, W2)|^ da;idw2 = 1, (A26) 


and given by [1 
4>{uJI,UJ2) = 


2\f2'Kc? i 


j5 hy v{uji)v{uj2) 

^ —r(/>p(a; 3 )(/)p(a; 4 ) 

y V{U}3)V(U}4) 

X S'(fc(wi),fc(w2),fc(w3),fc(w4)) 

X i5(a;i -I- a;2 — W3 — ^^ 4 ), 


daj3 / dw4 


(A27) 


where v{uj) = dui{k)/dk is the group velocity at frequency 
uj. The pump light is taken to be in a coherent state, 

(D (r,t)) = a J dk y^^/p(fc)DLr“(r)e““^('=)*-bc.c. 

(A28) 
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in the Heisenberg picture (compare (IA13|) '): in terms of 
the complex number a, the expectation value of the pho¬ 
ton number in the pump pulse is |a| if the function /p(fc) 
is normalized, 


J \fp{k)fdk = l. 


(A29) 


at w = wp and (/)p(wi + U !2 — strongly peaked at 
uji + UJ 2 — uj = ujp by the form (IHTl) adopted, we have 


(l){uJi,u;2) 


Z-Kiyf^a^h r WlijJ2 Wp 


4eo/3 Y v{uji)v{ui2) n(wp) 

X J(a;i, W 2 , wp, wp) J duj(f>p{uj) 


The spectral amplitude (j)p{ui)^ 


X 0p(cl’2 ^2 — ^)- 


(Bl) 


\(j)p{uj)\‘^ duj = 1, 


(A30) 


Using the spectral amplitude (IHTl) the integral over w can 
be performed using the identity 


is the function of frequency corresponding to fp{k), 
where from (IA29I) and ()A30ll we have |/p(fc)|^dfc = 
|0p(a;)|^daj. Finally, |/3|^ is the probability per pump 
pulse that a pair of photons is spontaneously created; it 
is determined by the expression (IA27|) for the BWF and 
the normalization condition (|A26p . 

Using (IA24|) in (IA16I) . substituting the result in (IA27I) 
and integrating over a;4, we find 

, , I ujiU!2 

4eo/3 \J v{uji)v{uj2) 

UJ{UJ1 + UJ2 - uj) , , . 

—r-7-^-r<(>p(w) 

v{UJpjv{UJi + UJ 2 — UJ) 

X (f>p[uJi -\- W 2 — io)J (wi, UJ 2 , w, UJi -h UJ 2 — Uj), 

(A31) 




where 

sin^ 

J (wi,a;2,W3,W4) . g 

sinf 

X k{^2), k{uj3),k{u4)), 

(A32) 

we use k{uj) again to indicate the wave number written as 
a function of frequency, and we now consider /r and x to 
depend on the frequencies {w„} through the dependence 
of the {fc„} on frequency; from (IA22I) . (IA25I) we now take 

fi =M({Wn}) 

= [k{uj3) + k{uj4) — k{uji) — fc(w2)] A 
-I- 9 [^(wa)] + 9 [fc(w4)] - 9 [fc(wi)] - 9 [k{uj2 )\, 

(A33) 


sinc(?/ — z) = — [ da:sinc(a; — y) sinc(x — z), (B2) 

Ti" J 


and we find 


(j){uJi,UJ2) 


dTTi'/2a^h I uJ]W2 wp 


4eo/3 

X J(a;i,a; 2 ,wp,a;p)sinc 


v{oJi)v{uj2) ^^(wp) 

2u;p — oji — UJ2 


Aui 


(B3) 


Within these approximations we determine |/3| from the 
normalization condition (1281) . For Aw small enough, in 
the expression that results we can put 


1 .of 2wp — Wi — UJ2 
sine ' 


ttAw 


Aw 


and thus the double integral in 
integral, giving 


,2 _ 27r 

' ^ AT 


37r-\/2a^fi 

w2 / 

4eo 

V^{ujp) J 


5(2wp — wi — W2), (B4) 
reduces to a single 

Wi(2wp — Wi) 

u(wi)u(2wp — wi) 


dwi 


X I J(wi, 2wp — wi, wp, wp)|^ . 


(B5) 


Recall that wi ranges over frequencies of the signal and 
idler. Because of the symmetry of the integrand, we can 
replace this by an integral that ranges over only the fre¬ 
quency range of the signal, if we simply multiply by 2. 
This yields (15^ . Within the same approximation the 
joint spectral density |(/)(wi,W2)|^ that follows from (IB3I1 
is immediately seen to be (1551) . 


X =x({w„}) 

= [fc(w3) -I- fc(w4) — fc(wi) — fc(w2)] Zi 

+ N{9 [fc(wi)] + 9 [fc(w2)]} + {N- l)/i/2. (A34) 

Appendix B: Long pulse limit details 

In this Appendix we construct the limit of the ex¬ 
pression (1551) for |/3|^ assuming a very long pump pulse. 
Taking ((55)) we note that for 4>p{uj) strongly peaked 


Appendix C: Connection to Fermi’s Golden Rule 

To connect the result for the photon generation rate in 
the long pulse limit (I34|) to the simple result of Fermi’s 
Golden Rule m we need to relate the am (1551) associ¬ 
ated with the individual rings to the amplitude a char¬ 
acterizing the photon number in the pump pulse appear¬ 
ing in (IA281) . That is, we must identify how the inci¬ 
dent pump pulse effectively “loads” each ring with pump 
light. And to do that we must first relate the mode held 
Diso(r, m) of m characterizing an isolated ring m to 
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the mode field D™®(r;m) characterizing ring m in the 
asymptotic-in state (lAlOl). F or the convention adopted in 
m) satisfies la the normalization condition 


Drso(i', m) ■ Diso(r,m) 
eoe(r) 


dr = 1, 


(Cl) 


where the integral is over the m*** ring. In the approxi¬ 
mation (03 of isolated rings the pump frequency in the 
ring is wp, and letting kp identify fc(a;p) we can write the 
field associated with ring m in the asymptotic-in state at 
ujp (see (lAllM as 


D™s(r;m) [T(fcp)]’"-^ D™«(r-z^z) 

[T(fcp)]’” ^5'(wp) 




X>fcp(r - ZmZ) 


,ikp(^ 


(C2) 


Dr(r; 


where we have used ((Xill and the field enhancement fac¬ 
tor 5'(u;p) is given by (IA5I) . 

We want to identify the proportionality between 
m) and Diso(r, m), putting D^‘p®(r;m) = 

CDi,o(r, m), where (7 is a constant to be determined, 
r • 1 * 

We integrate D^'"®(r;m) • D™®(r; m)/(eo£(r)) over 

the rn}^ ring; that integral will involve an integral over 
the length of the ring, which will give a factor I = 2'kR 
with I the circumference of the ring and R its radius, as 
well as an integral over the plane perpendicular to the di¬ 
rection of propagation of the light in the ring. Recalling 
that we have taken Rkp (r) to describe the same physical 
mode profile as dtp (rj_) (see text after (IA4I) 1. we can use 
the normalization condition [s^ of dfc(rj^), 

f dJ(rj-) • dfc(rj_) 


to 


eoe(i') 

the integral 


-dxdy = 1, 


(C3) 


of 


D^,7(r;m) 


complete 

(r; m)/(eoe(r)) over the ring. Using (ICll) we 
can then identify the constant C, within a phase, and 
write 

D™®(r;m) = [r(fcp)]’"-^ 5 (c^p) Vi? Di,o(r,m). 

(C4) 

The phase chosen is (except for perhaps a global phase 
factor) the natural one, since it is natural to take 
Diso(i‘,m) = Diso)!" — ZmZ), where Diso(r) would be the 
mode field of an isolated ring with its closest point to the 
bus waveguide at z = 0; that is, the Diso(i‘,m) are es¬ 
sentially all the same, except for the position of the ring, 
and so the phase factors associated with the propagation 
of the asymptotic-in state will reside in D^‘p®(r;m). 

We can now identify the a^. of (l2^ that would be 
associated with the pump field in each ring for excitation 
by a long pulse characterized by a. From (IA28I) we have 


(D (r,t)) =c 


4>p{u!) 


dk hw duj 

dw V 2 V dfc 

-I- C.C., 


-l-jasy-m/ N 


(C5) 


where we have converted to an integral over frequency 
and used |/p(fc)|^dfc = \(j)p{uj)\^ duj (see text after 

(IA30() 1. For a <))p(w) strongly peaked at ujp we can write 
DLfp”(r), and put dw/dfc « w(a;p) and 
Lj Ri wp in the prefactors of (j)p(uj) to write 






, dw 


+ c.c. 




Trhujp 


D 


asy-in 


(r)- 




'j{ujp) 

AT AT 

for-— <t< ——, 

2 2 

=0 otherwise, 


+ c.c. 


(C6) 


where in the last two lines we have specialized to the 
pump pulse of (1501) . Now using the expression (lAlOl) for 
DL^^““(r), and (IC4I) for the field in the ring, we find 


(D (r,t))ri„ 




X Diso(r, m)e -I- c.c. 

AT AT 

for- <t < - 

2 2 

iO otherwise. 


(C7) 


Since from the expression (1141) for a pump field loaded 
into the rings we would have 


(D (r,t))^ 






;r,m) 


^ — iujpt 


+ C.C., 

(C8) 

for the time AT that we were making a calculation using 
Fermi’s Golden Rule, we can identify the relation 


I 


v{u!p) (AT) 


g*fep.™[T(fcp)]’"-^5(c^p), (C9) 


where again I = 2ttR is the ring length and R the ring 
radius. 

To now make the connection of the result (IMl) with the 
calculation using Fermi’s Golden Rule, we look at 


a^J(wi, u;2,t^p,wp) =a^ ^ / drF*'^j'(r) 

m 

* [of; 

X (cw) 


),asy-out 


(3) 

y.asy-out 
Z,asy-in 


(l‘)^Rfe 


where we have used (IA32I) together with (IA24p and the 
definition (IA17|) : here fci = k{u!i) and ^2 = k{uj 2 ), the 
sum is over all the rings, and the integration over volume 
Vm over is over ring m. Within ring m we have D“V(™(r) 
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equal to D^™®(r;m), and using the result (IC4I) we have 


J(a;i,a;2,w,a;i +W 2 — uj) as 


J(a;i, a;2, wp, wp) 

= [T(fcp)] 


2(m-l) ^2 


;?^(c.p)i7/ drr;^3;'(r) 


m 

v{wp){AT) 


271 


■J2alF(^Hu;uco2), 


(Cll) 


[fc(wi), k{u}2), k{uj), fe(wi + a;2 - w)] 

_ / 2 y A/2 A/2 

VI - (T j (‘‘^I - ^s) - i(A/2) (a;2 - wi) - i(A/2) 

^ A/2 _A/2_ 

(w — wp) — i(A/2) (w — wi — a;2 — wp) — i(A/2) 

X y'^r(y(r)[P*(r)P^(r)]*P'=(r)P'(r), (Dl) 

showing the expected resonance behavior. From the gen¬ 
eral expression (l36t for J(a;i,a;2, W3,a;4) we now seek 


where in the third line we have recalled the definition 
(HH) of W2), and in the fourth line we have used 

the relation (IC9I) between am and a. Using (jCll|) in 
(l34)l we find 


\£l 

AT 


^ dirh^uip 
8^2 


dw 


w(2a;p — oj) 
v{oj)v{2ojp — ui) 


^a^uW(u;,2^P-a;) 


(C12) 


To compare this with the result (EU from Fermi’s 
Golden Rule, in that latter result we pass to an infi¬ 
nite normalization length £; the distance between the 
rings and the number of rings of course remains fixed; 
then ks and ki become continuous variables, j,. —>■ 

{C/2tt)^ f dfcidfc2, and integrating over all possible final 
states we have 


J(a;i, W2, w, wi -I- W2 — w) 

• uAT 

sm ^ 

= . I f'‘\k{wi),k{u}2),k{u}),k{u;i + u}2 - w)), 

(D2) 

where the general expression (IA33I) for ^ reduces, in the 
limit (1441) of vanishing group velocity dispersion to (recall 

(HTt and (1501) 1 

^ = 9 [A:(a;)] -I- 6 [fc(wi + UJ2 — ^)] — 9 [A:(a;i)] — 9 [k{uj2)] , 

(D3) 

and the general expression (|A34I) for x reduces in that 
same limit to 

X = N{9 [fc(a;i)] + 9 [k{u;2)]} + - l)/2. (D4) 

Here, for simplicity, we have put zi = 0 and, in both of 
these expressions (see m). 


m 

t 


QTT^^O/p 

8eg 


dfcidfc 2 UI1L1J2 


^a^FM(a;i,i^2) 


X S{huJi + hi02 — 2Hu!p). 


(C13) 


Changing the integrals to those over frequency, and using 
dk/du! = l/v{u!), where v(uj) is the group velocity, we 
can integrate over da;2 with the use of the Dirac delta 
function and find 


m 

t 


9Trh?Ljjp 

8^2 


dw 


w(2a;p — w) 
u(a;)u(2a;p — w) 


Y,al,F^^Hu;,2u;p-u;) 


(CM) 


in agreement with (IC12I) . 


9{k) = TT + kl + 2 tan ^ ^. (D5) 

\ 1 — cr cos kl) 

In the limit of weak coupling (HOI) and for small deviations 
in k (w) from a resonance, we simplify this as 

9(k) Ki TT + kl + 2 tan“^ . (D6) 

Then, again neglecting group velocity dispersion and 
Taylor expanding fc (w) to first order (HOI) , as well as re¬ 
calling the resonance condition (HD) , we write 

9 [k{uj)] Ki{2M+l)TT+ ~ 

+ 2tan-i(^^t^). (D7) 


Appendix D: The coherence number 

With frequency ranges restricted according to (ISD) and 
|u| < i5 (recall (|52l) l. here we calculate a useful approx¬ 
imation for the number of rings that will behave coher¬ 
ently. We first rewrite the single ring component (H51) of 


For the first two terms of (imi) . the pump terms, i5 <C A 
guarantees that the argument of tan“^ is small and we 
have 

9 [k (w)] + 9[k (uji +UJ2 — w)] 

f«2^(2P+l) + - + ^. (D 8 ) 

Va A 
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For the second two terms of (imi) . the generated 
terms, we cannot simplify the tan“^ functions the same 
way, since from the resonant behavior in (IDll) we 
must consider wi — ws and 102 — ranging over A. 
However, using the identity tan“^ (x) + tan“^ (y) = 
tan“^ [(x + y) / (1 — xy)] we find 

6[k (wi)] +6[k {UJ 2 )] 

- 27r(5' + /+ 1) + — +2tan~^ ^ ^ , (D9) 

where we have dehned 

rj = [lui - ujs) - iuj 2 - LOi). (DIO) 

The denominator in the tan“^ function of (IDOp is always 
positive and of order A or larger, since |u| is of order 
S and 6 A, and so the argument of that function is 
always much less than unity; keeping the first term in the 
expansion we write 


last term in (IDllI) and approximate A^/(A^ + ?7^) fa 1/2. 
Then taking a typical value of |u| to be about 5, we char¬ 
acterize the typical size of the expressions in (ID8I) and 
(IDllI) to be 

9 [k (w)] + 9[k {oji + 0 J 2 — w)] 

« 2^ (2P + 1) + <5, (D12) 

and 

9[k {LOi)] + 9[k (a;2)] 

fa 2 tt{S + I + 1) + (J-+ S (D13) 

respectively. From (ID3I) . this gives a typical size for y of 
about (1551) . We also note that a typical size of the phase 
X appearing in (ID2p for J{uji,uj 2 ,uj,uji -|-W2 — oj), given 
by (ID4I) . is found from (ID13I) and (1^ as 


9 [k (wi)] 9 [k (^ 2 )] ~27r(5' -I- / -I- 1) H - 

Vg 

4u A^ 

A A^ -f _ y2 


(Dll) 


For a pump pulse with \cj)p{uj)\^ peaked at the ring res¬ 
onance frequency wp, energy conservation and the reso¬ 
nance structure of j Im]) will allow for a range of oJi 
and UJ 2 with ry varying over a range of about —A/2 to 
A/2. For most of this range we will have 77 ^ ^ so as 
a rough estimate we neglect in the denominator of the 


X=2,iV(S + /+l)+(lv(F + |)+^)i 

f^27TN{S + I+l)+(^^y (D14) 

where in the second line we have used our assumption 
that ^/A <C 1 and have taken the finesse of the rings to 
be very high, so that the free spectral range (1551) is much 
larger than the width A of the ring resonances. Under 
these conditions the phase x will vary significantly over 
the BWF only as N becomes a significant fraction of N^oh 

dSSl). 
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